Tunable Charge and Spin Seebeck Effects in Magnetic Molecular Junctions 
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We study the charge and spin Seebeck effects in a spin-1 molecular junction as a function of 
temperature, applied magnetic field, and magnetic anisotropy (D) using Wilson's numerical renor- 
malization group. A hard- axis magnetic anisotropy produces a large enhancement of the charge 
Seebeck coefficient S c (~ fcs/|e|) whose value only depends on the residual interaction between 
quasiparticles in the low temperature Fermi- liquid regime. In the underscreened spin-1 Kondo 
regime, the high sensitivity of the system to magnetic fields makes it possible to observe a sizable 
value for the spin Seebeck coefficient even for magnetic fields much smaller than the Kondo temper- 
ature. Similar effects can be obtained in Ceo junctions where the control parameter, instead of D, 
is the gap between a singlet and a triplet molecular state. 

PACS numbers: 72.15.Qm, 72.15.Jf, 73.63.-b, 75.76.+j 



The Seebeck effect refers to the generation of a charge 
current (or a voltage drop) by a temperature gradient 
applied across a metal The spin- Seebeck effect [2|-[4| , 
concerns the thermal generation of pure spin currents. 
Applications of the Seebeck and spin-Seebeck effects at 
the nano-scale, with potential impact on a variety of new 
technologies, could profit from the scalability and tun- 
ability properties of nano-devices but still require a bet- 
ter understanding of thermopower effects in nanostruc- 
tures and of the effect of strong electron-electron corre- 
lations on them. The recent experimental observation of 
the Seebeck effect in different nano-structures, in partic- 
ular in molecular junctions [5J and quantum dots (QDs) 
0, opened new routes to study these phenomena. Here 
we show that junctions with spin-1 molecules, like the 
Co(tpy-SH) 2 complex [7] or the Ceo buckyballs [8], give 
rise to a large and controllable enhancement of the See- 
beck and spin Seebeck effects at low temperatures and 
low magnetic fields. We also show that thermoelectric 
experiments in these systems give access to valuable in- 
formation on the residual interaction between quasipar- 
ticles in the low temperature Fermi-liquid regime |9j, in- 
formation that cannot be obtained by conventional tech- 
niques. This offers the opportunity to experimentally 
study these effects on the smallest possible lengths scales 
looking for the basic mechanisms of thermomagnetic ef- 
fects in strongly correlated systems flol-[l3| — in particular 
in the regime where energy transfer is governed by spin 
fluctuations. 

We start by describing the case of junctions with the 
organometallic molecule Co(tpy-SH)2 [see Fig. 1]. In this 
complex, the Co ion is located in an approximately octa- 
hedral environment and is found to be in the Co 1+ 3d 8 
configuration with a total spin 5=1. When the molecule 
is mechanically stretched, the local environment of the 
Co ion is distorted and a magnetic anisotropy, with an 
easy plane, emerges due to the spin-orbit interaction. As 
the molecular orbitals reproduce the magnetic structure 
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FIG. 1. Schematic representations of the Co(tpy-SH)2 and 
Ceo molecular junctions. In the former, a distortion of the 
crystal field and/or an external magnetic field induces the 
splitting of the 5=1 state. In the Ceo case there is a singlet- 
triplet gap given by the exchange constant J. 



of the metallic ion, the molecular junction is modeled by 
the Hamiltonian H = Hm + ^iv + T~Li where 



H M = ^2 [^t%+ £ (^t + %)] 

£=a,b 

-gii B BS z + DS 2 z , 



JS a • Sb 



(1) 



describes two degenerate orbitals (a and b) of the 
molecule with level energy e and Coulomb repulsion U, 
coupled through a spin exchange interaction J. Here 
S = S a + Sb is the spin operator, Si = J2a*> d la^^ ,d £a 
where & aa > the Pauli vector, and fi£ a = d\ a d t(J is the 
number operator. The molecular spin is coupled to an 
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external magnetic field £>, and D is the strength of the 
magnetic anisotropy. In this case, the orbitals describe 
the Cobalt e g levels and J < is a Hund rule ferro- 
magnetic coupling — in the Cqo case, J > is antiferro- 
magnetic and D = 0. The hybridisation between the 
molecular orbitals and the leads' states is described by 

= Ekaa Vka (4a^kaa + C Loc d aa) > where a = L,R 

stands for the left and right electrodes, respectively. Here 
we have assumed that a single conduction channel is ac- 
tive in order to describe the underscreened Kondo effect 
observed in Ref. [7j, and we have considered, for sim- 
plicity, that only one of the molecular orbitals is coupled 
to the electrodes. The Hamiltonian of the metallic leads 

is given by U\ = E a k* £ ka C La C k*a' We focus on the 

regime e < 0, and U ^> \s\ where only the single and 
double occupied states of the molecule are relevant for 
the low-temperature physics. 

The presence of a voltage (V) and temperature drop 
(AT) across the junction generates a charge current that 
in the linear response regime is given by 



GV + GS C AT 



(2) 



where the conductance G and charge Seebeck coefficient 
S c are given by 



G 



S r = - 



k f 



el k B TX 



(3) 



with l n = £ CT 1° and 2£ = " n (- ^) %(u>) dw . 
Here the quantity T a (oj) describes the tunneling of spin- 
a electrons across the junction and is given by T a (oj) = 
ri+rfl ^^) where /v(co>) is the spin dependent spectral 
density of the molecular state. In the above expression 
Tq; = npaVa is the contribution to the width of the 
molecular energy levels introduced by the coupling with 
the a lead, p a is the electronic density of states per spin 
of the electrodes at the Fermi level and f(uS) is the Fermi 
function. In what follows we assume p a to be constant, 
set the Fermi energy to zero (ep = 0) and choose the 
bandwidth of the leads W as the unit of energy. 

The molecular spectral density is evaluated using 
the numerical renormalization group (NRG), a non- 
perturbative technique known to give excellent results 
for the transport integrals [14|, LL5j. The total conduc- 
tance and the charge Seebeck coefficient S c are shown in 
Fig. [2] for different values of the anisotropy parameter D 
and B — 0. For the isotropic case (D = 0) the molecular 
spin is partially screened by the conduction electrons be- 
low a Kondo temperature Tk (the under-screened Kondo 
effect) with Tk ranging between IK and 200K depend- 
ing on the device [7]. The ground state of the system is 
a singular Fermi-liquid with a free (unscreened) spin \ 



Ifil El- The development of Kondo correlations is asso- 
ciated with a monotonic increase in the conductance as 
T is lowered (top panels Fig. [2}. With a spin anisotropy 
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FIG. 2. Left panels: color maps of the conductance G (top) 
and charge Seebeck coefficient S c (bottom) as a function of 
the temperature and the anisotropy constant D. For finite 
D, the conductance goes down to very small value s for T < 
T^ while the S c reaches its maximum, (fcs/|e|) ir/^/3(l + 7), 
for T < Tk- Right panels: temperature dependence of G 
and S c for different values of D, indicated in the color maps 
with dashed lines. The conductance is given in units of go = 

(2e 2 /h)4r L r R /(r L + r R ) 2 . 



< D < ksTx , as the temperature is lowered the con- 
ductance first increases following a universal behavior, 
goes through a maximum and decreases at a characteris- 
tic temperature T£ = T K exp(-2 ^k B T K /D) 0. This 
new temperature scale can be identified with a second 
stage Kondo effect, induced by the magnetic anisotropy, 
in which the remaining spin \ is screened. In this regime 
the ground state of the system is a Fermi liquid. 

Within this scenario, that has been confirmed by ex- 
periments and theory, the Seebeck coefficient shows novel 
features (bottom panels Fig. [2]). For D = it develops 
a large peak at very high temperatures (unphysical for 
magnetic molecules) and changes sign at a lower temper- 
ature. This behavior is similar to what is obtained for 
spin- ^ quantum dots 

nnm. F ° r the stretcned (d > 0) 

molecule [19], a second large peak develops at T < T£. 
This low-T peak is controlled by spin fluctuations and it 
is strongly dependent on the spin anisotropy. For small 
anisotropy (D < ksTx), the two characteristic energy 
scales T£ and Tk are far apart and S c is a function of 
T/T K at low temperature. For larger anisotropics, the 
scaling breaks down and the maximum value of S c de- 
creases. In the scaling region, and for T/T K < 1, the 
system is a Fermi-liquid and S c can be described by a sim- 
ple analytical expression. In Fermi-liquids, the residual 
interaction between quasiparticles leads to a self-energy 



with an imaginary part proportional to (uj 2 + 7r 2 T 2 ) 
[I, [13]. Using this result to estimate the spectral den- 
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FIG. 3. Spin Seebeck coefficient S s and spin current I s for 
an isotropic molecule (D = 0). Top panel: S s = St — S^ 
as a function of the magnetic field for different temperatures. 
Note that S s is significantly large even for g\iBB <C ksTx- 
The maximum of S s occurs for g\iuB ~ ksT. Bottom panel: 
Spin current as a function of gfiBB in the zero charge current 
condition (I c = 0). 



sity p a (w) we obtain [2l|, 

k B 27r 2 (k B T/uj ) 



S c = 



^(l + 7 )(fe7>o) 2 



(4) 



Here ooq = fc^Tj^ sin(7rn/2), n is the total occupation of 
the molecular orbitals and 7 > measures the strength 
of the quasiparticle interaction. Notably, the maximum 
of the Seebeck coefficient 5™ ax = (k B / \e\)n/^3(l + 7) 
depends only on universal constants and the parameter 7 
that reduces the maximum value of the thermopower — 
the coupling to the leads, D and other parameters that 
are both experimentally hard to determine and sample 
dependent, do not affect the value of 5™ ax . The tem- 
perature behavior of the Seebeck coefficient is one of the 
central results of this work. Its measurement provides 
a direct access to the parameter 7 and can be used to 
quantitatively test Nozieres local Fermi liquid ideas [9j. 

The fitting of the numerical results with Eq. (j4]) is 
excellent (see [21[) and gives 7 ~ 2.8 for the parameters 
used in our simulation. Therefore, in a junction with 
Tk — 100K, the Seebeck coefficient could be tuned to 
have a peak of the order of k B /\e\ ~ 100 fiV/K at a tem- 
perature in the range of 0-50K. 

Let us now focus on the effect of an external magnetic 
field. There are now two control parameters, the mag- 
netic field B and the anisotropy D, leading to different 
magneto-thermal regimes for the junction. Defining the 



contribution to the charge current of the spin-a electrons 
as la = G(jV -\-GaScrAT, where G a is the a-contribution 
to the total conductance [cf. Eq. (j3j)], the spin dependent 
Seebeck coefficient results, 



|e| k B TX% " 



(5) 



If the condition I c = is set, a voltage drop across the 
junction will develop as before but in addition, the tem- 
perature difference AT will also generate a pure spin cur- 
rent i s = (/t - h)/\ e \ & iven b y 



(fit - s±) 



AT 



(6) 



with 1Z = G a 1 . In this units I s gives the number of 
spins per second flowing though the junction. Note that 
with these definitions S c = T,a( G °/ G ) S v For D = °> 
the under-screened Kondo regime, the system behaves as 
a singular Fermi-liquid with a divergent spin susceptibil- 
ity as T — > [18] and p(uj) shows a Kondo peak with a 
singular behavior at the Fermi energy . As the ground 
state of the system contains a free —spin, the presence of 
a magnetic field is a relevant perturbation in the renor- 
malization group sense. Consequently, even a very small 
magnetic field is able to spin polarize the molecule and 
split the Kondo peak in p(uj) (see [ll[). The magnetic 
field dependence of S s = — S± and I s is shown on Fig. 
03 This spin polarization of the molecule results in the ap- 
pearance of a pure spin current even for g\i B B <C k B Tx 
(but g\i B B ^ k B T), in clear contrast with the S = \ 
case (Anderson model) where g\i B B ~ k B T^ is required 
in order to observe a sizable spin polarizing effect [12j . 

Similar results are obtained for the stretched molecule 
only if g\i B B > k B T^. For g\i B B < k B T^ the coefficient 
S s is small and a large charge Seebeck peak is recov- 
ered (see Fig 3) This behavior shows that with some spin 
anisotropy (D ^ 0) the voltage drop V at the junction is 
very sensitive to the external field and can evolve from a 
large value (of the order of k B AT/\e\) for B = to zero 
for g\i B B > k B T^ [22[. The low temperature behavior 
of the spin dependent Seebeck coefficients S a can also be 
described using a local Fermi liquid theory. We obtain 
expressions similar to Eq. ((4]) where now ujq is replaced 
by oJoa = k B T^ sin(7rn a -), with n a = n/2 — am/ 2 the to- 
tal number of spin-a electrons in the molecular orbitals. 
As n ~ 2, a small magnetization m can change the sign 
of one of the energies uo a and consequently the sign of 
one of the spin dependent Seebeck coefficients S a leading 
to a large spin current. 

Ceo junctions present similar thermopower properties. 
A Ceo molecule embedded in a junction may have a sin- 
glet 5 = state and a triplet 5=1 state separated 
by a small singlet-triplet gap J that can be controlled 
by gate voltages |8]. When the singlet state lies just 
below the triplet state the system shows a two stage 
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FIG. 4. Left column panels: color maps of the Seebeck coefficient (S c ), the spin current (J s ) and the spin dependent Seebeck 
coefficients (S a ) as a function of the temperature and the anisotropy constant D. In all cases g\±BB — 0.01/cbTk Note that 
S c is strongly reduced for D/Tk < 0.5 that corresponds to g\iBB < ksT^ and that I s reaches its maximum value in the 
parameters region where S± changes sign. Notice that, as D decreases, the S\, peak moves down to very low T before changing 
its sign. Right column panels: Temperature and anisotropy dependence of S c and I s . Note that S c is suppressed in the region 
where I s is different from zero. 



Kondo effect with a Kondo temperature Tk, of the or- 
der of 4K in [23| , and a second stage temperature given 



by Tj£ oc T K exp(-k B T K /J) [23|, |24|. While the global 
behavior of the conductance is similar to that of the 
anisotropic spin-1 case, its specific high T-dependence 
is different. By contrast, the temperature dependence of 
S c for T < ksTx is the same as the one obtained for 
the Co-complex and it is well described by Eq (j4]) with a 
particular value for 7 that depends on the residual inter- 
action in the Cqo case. If the triplet is the ground state, 
and for J > Tk, the system is described by the isotropic 
under-screened S = 1 Kondo model and the physical be- 
havior corresponds to the D = case discussed above. 
It is worth mentioning that other systems, like an artifi- 
cial molecule build with a T-shape double Quantum Dot 
(2QD) may also show some similar effects. In these two 
cases, however, the Kondo temperatures are usually very 
small making it more difficult to experimentally observe 
them. In this sense, the stretched Co complex presents 
advantages over these systems that make it an excellent 
candidate to observe high charge and spin Seebeck effects 
at low or moderate fields and temperatures. 

In summary we have presented a study of the charge 
and spin Seebeck effects for real molecular junctions that 
combines the 5 = 1 underscreened Kondo effect together 
with a controllable spin anisotropy Dora single-triplet 
gap J. We found that these systems presents a large 
charge Seebeck coefficient 5 C , of the order of fe^ a ^ a 
temperature that can be controlled either by stretching or 
gating the molecules. Quite remarkably, the maximum 
value of S c in the low T regime, (&B/|e|)7r/y3(l + 7), 
depends only on the Fermi liquid properties of the sys- 
tem (characterized by 7). It is important to emphasize 



that this low T behavior is fully controlled by spin fluc- 
tuations. In the presence of a magnetic field these sys- 
tems are able to filter the spin flow in a regime where 
Kondo correlations are fully developed and present mod- 
erate values of S s and I s even for magnetic fields such 
that the Zeeman energy is much smaller than the Kondo 
energy scale. 
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